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Abstract: The problem of feedback control of quantum systems by means of weak measure- 
ments is investigated in detail. When weak measurements are made on a set of identical 
quantum systems, the single-system density matrix can be determined to a high degree of 
accuracy while affecting each system only slightly. If this information is fed back into the 
systems by coherent operations, the single-system density matrix can be made to undergo 
an arbitrary nonlinear dynamics, including for example a dynamics governed by a non- 
linear Schrodinger equation. We investigate the implications of such nonlinear quantum 
dynamics for various problems in quantum control and quantum information theory. The 
nonlinear dynamics induced by weak quantum feedback could be used to create a novel 
form of quantum chaos in which the time evolution of the single-system wave function 
depends sensitively on initial conditions. 

The conventional theory of quantum feedback control assumes the use of strong or 
projective measurements to acquire information about the quantum system under control 
(1-10). Such measurements typically disturb the quantum system, destroying quantum 
coherence and giving a stochastic character to quantum feedback control. But strong 
measurements are not the only tool available for acquiring information about quantum 
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systems (11-12). In nuclear magnetic resonance (NMR), for example, one makes collec- 
tive measurements on a set of effectively identical systems: by monitoring the induction 
field produced by a large number of precessing spins, one can obtain the average value 
of their magnetization along a given axis while only slightly perturbing the states of the 
individual spins (13). We will call such measurements 'weak measurements' since they 
only weakly perturb the individual systems in the set. (Such weak measurements on large 
sets of identical systems should not be confused with the single-system weak measurements 
debated in (14-16)). The information acquired by weak measurement can then be fed back 
to the spins, for example to suppress super-radiant decay (17-18). NMR is not the only 
system in which weak measurement is possible: one can perform weak measurements on 
essentially any set of quantum systems that can be coupled weakly to an external ap- 
paratus. This paper provides a general theory of quantum feedback control using weak 
measurements. Since weak measurements allow the accurate determination of the com- 
plete single-system density matrix of each member of a set of identical quantum systems, 
while effecting each system in the set arbitrarily weakly, quantum feedback by weak mea- 
surement will be shown to be capable of accomplishing tasks that are not possible using 
conventional, strong measurements. A model of quantum feedback using weak measure- 
ments is given and applications are proposed. In addition to NMR, quantum feedback by 
weak measurements could be used in quantum optics and atomic and molecular systems to 
effect arbitrary nonlinear Schrodinger equations, to create solitons and Schrodinger cats, 
to perform quantum computations, and to institute novel forms of quantum chaos. 

Quantum feedback via weak measurement represents a novel paradigm for coherent 
control of quantum systems. It allows the performance of operations that are impossible 
in the normal, strong measurement paradigm for quantum control. For example, suppose 
that each of the systems in the set is in the same unknown pure state. Then feedback 
with weak measurement can be used to drive them to any desired pure state reversibly, 
while preserving quantum coherence. This contrasts markedly to quantum feedback using 
strong measurements, where a system in an unknown quantum state can be driven to any 
desired quantum state, but only at the cost of disturbing the system's state irreversibly 
and stochastically, destroying quantum coherence in the process. 

The general picture of quantum feedback control using weak measurements is as fol- 
lows. Suppose that we have an set of N identical non-interacting quantum systems, each 
with density matrix p. (Of course, no set of quantum systems is perfectly non-interacting. 
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but in many situations — e.g., liquid-state NMR, quantum optics — the non-interacting 
approximation holds to a high degree of accuracy.) Assume that the system is coherently 
open-loop controllable, so that we can perform arbitrary unitary transformations U on the 
system (necessary and sufl&cient conditions for open-loop coherent control of quantum sys- 
tems are well known (1-6)). Now assume that we are able to make a sequence of collective 
weak measurements on these systems that allow us to determine the single-system reduced 
density matrix p to some degree of accuracy 6, while disturbing this density matrix by an 
amount e. As will be seen below, both 6 and e can go to zero in the limit that the number 
of systems goes to infinity. If the systems are individual nuclear spins, for example, the 
single-spin density matrix can be determined by measuring the induction signal produced 
about two different axes: this allows one to determine the expectation of the magnetiza- 
tion along the x, y, and z axes, which is in turn sufficient information to determine the 
single-spin density matrix. Now feed that information back into the set by applying to each 
system a unitary transformation Us{p), where Us is some potentially nonlinear function of 
p, and the subscript S indicates that Us discriminates between different p to an accuracy 
5. The time evolution of the system with feedback by weak measurement is accordingly 
given by 



where Ap is the perturbation to the single-system reduced density matrix induced by the 
weak measurement, with ||Ap|| < e for a suitable norm || || such as the sup norm. As will 
be shown below, in the limit N ^ oo, the collective measurement can be performed in 
such a way that both S and e ^ 0, and the time evolution of the single-system density 
matrix is governed by the equation 



The remainder of this paper will be devoted to exploring the implications of equations 
(la) and (lb). These equations have a variety of interesting features. The first, perhaps 
most obvious, is that they can be nonlinear as a function of p: if p = api +(3p2, it need not 
be the case that U{p)pU\p) = aU{px)p\U^{p\)+l3U{p2)p2U'^{p2)- (It is important to note 
that although the single-system reduced density matrix undergoes a nonlinear evolution, 
the density matrix for the set of systems taken collectively undergoes a conventional linear 
time evolution: no laws of quantum mechanics are broken in constructing this nonlinearity.) 



p' = C/,(p)(p+Ap)C/](p) 



(la) 



p' = U{p)pU\p) 



ilh) 



3 



If the weak measurement is made continuously in time, then in the hmit N — > oo, 5 — > 0, 
e — > 0, feedback causes the single-system density matrix to obey a nonlinear Schrodinger 
equation 

dp/dt=-i[H{p),p] (2) 

where H{p) is the Hamiltonian corresponding to U{p). Such nonlinearities in the case of 
sets of nuclear spins are well-known: for example, if each nuclear spins in the set interacts 
with the mean field generated by the spins taken together, then the single-spin density 
matrix obeys a nonlinear Bloch equation (19). Feedback by weak measurement allows 
one to impose an arbitrary nonlinear Hamiltonian dynamics on the single-system density 
matrices: if in the open-loop case, without feedback, one can apply any conventional 
linear time evolution, then in the closed loop case, with feedback of the results of weak 
measurements, one can apply any desired nonlinear dynamics that preserves the eigenvalues 
of the density matrix. That is, one can take p f{p), where f{p) has the same eigenvalues 
as p. If one can apply open-system operations (20-21) as well as closed-system, unitary 
transformations, then one can alter the eigenvalues of the density matrix as well to take 
p — > g{p), where g{p) can be an arbitrary density matrix. 

Now let us look more closely at the dynamics of the weak measurement process, in 
order to determine how accurately the single-system density matrix can be measured and 
at what cost. There are two measures of the cost of weak measurement: first, the size N 
of the set required to attain a given accuracy S, and second, the amount e by which the 
individual systems are perturbed by the weak interaction with the measuring apparatus. 
Here we construct a specific model of weak measurement applicable to a wide range of 
physical systems. 

The general picture of measurement on N identical systems is as follows. The density 
matrix for the systems is ptot = p ® p ® ■ ■ ■ ® p- A positive operator valued measure 
(POVM) on this system corresponds to a set of operators {A^^} such that ^jii^M* ~ 
where I is the identity operator; the measurement corresponding to the POVM gives the 
result p, with probability = tr A^iPtot^jti^ in which case the system is left in the 
state ptot/u = i^/Pfi) Si ^M»/^tot^|ti and the density matrix for the fth subsystem goes to 
Pefj, = tr^/^^ptotM- 

We will define a weak measurement to be one that leaves the single system density ma- 
trices unchanged to within a small accuracy e: X^^P^uHp — pe^W < e- For example, a useful 
POVM is the set of Gaussian quasi-projections: A^, = (l/(27r)V^AV2) e-^""'^)'/'^^' |a)(a|f^a| 
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where the normahzation is chosen so that trAj^^An = 1 and Aj^Afj^d/j, = I (here there is 
no need for the auxihary index i). If we write the single-system density matrix in the a basis 
^ X^aa' ^aa'lci^ici'l, then the measurement corresponding to the determines the value 
of a = tipA to an accuracy A, where A = J a\a){a\da. In addition, the measurement has 
the effect of reducing the off-diagonal terms of p by an factor ^-ia-a')^ /2A'^ ^ corresponding 
to a perturbation of size e fti AvlV2A2, where A A = ^tvpA'^ - a^. If A » AA the mea- 
surement perturbs the system only weakly. Of course, the more weakly the measurement 
perturbs the system, the less information it acquires. By making a weak measurement on 
all the systems in the set simultaneously, however, one can obtain very precise information 
about the single-system density matrix while perturbing it only slightly. Consider the 
A'"-system POVM given by 

J —oo 

If a collective measurement corresponding to this POVM is performed on the systems in the 
set, one obtains the value of d to an accuracy A'^/N'^ + AA^ /N^ while still perturbing 
the single-system density matrix by the amount e ~ Ayl^/2A^. It can be clearly seen that 
in the limit A?" — > oo we can take A oc y/N, giving an arbitrarily accurate determination 
of d together with an arbitrarily small perturbation of the single-system density matrix. 
After the measurement, the over-all density matrix is in the form p (8) . . . <8) p + 0(e), so 
that the assumption of no correlation between the systems is only true to order e. In the 
limit A?" — > oo, e — > 0, however, the no-correlation assumption still holds. 

Now we construct a model of how such a weak measurement might be performed. 
Our model is analogous to weak measurements in NMR, in which each system in the 
set is coupled weakly to the electromagnetic field in the measurement coil. Couple each 
system to the measurement apparatus via a single continuous quantum variable ('pointer 
position'), described by an operator Q — J q\q){q\dq, via a Hamiltonian coupling jAP, 
where P is the momentum corresponding to Q; [P, Q] = i. This gives a dynamics for the 
system and pointer: \a)\q) |a)|g-|- 07^) over time t. Now suppose that all N systems 
are coupled symmetrically to the pointer by an interaction {Ai -|- . . . + An)P. If the 
systems are all originally in the state = p as above, and the pointer is originally in the 
state lip) = J ip{q)\q)dq, then the interaction between the systems and the measurement 
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apparatus gives 

pi (g) . . . (g) Pat ® PSM{t) 



^ / dqdq''ip{q)'ip{q')pa^a[ ■ ■ ■ Pa^a'^ 



(3) 



OlOj...OjvO^ 

\ai){a\ \ (8) . . . (g) |aAr)(a'Arl ®\q+{a\ + + aN)lt){q' + (a'l + . . . + a'jv)7t|. 



One can then find the state of the apparatus at time t by taking 
pM{t) =tr5 pSM{t) 

= Yl J dqdq'^p{q)^p{q')paiai---PaNaN (4) 

Ol ...OiV 

\q+{ai + ... + aN)lt){q' + (oi + . . . + aAr)7^|- 

That is, after the measurement the pointer registers the sum of N independent samples 
of A, where each result a occurs with probability pa = Paa- The sum is registered to an 
accuracy AQ = ■\/ {'4>\Q'^\'4') ~ (V'lQlV')^- measures the initial spread of the pointer 
wave function Accordingly, after the coupling of the pointer to the systems, the 

pointer registers the result {A) = irpA ± AA^/T/eN^~+T/N , where e = {'jt/AQ)'^ wiU 
be seen to be a measure of the degree of perturbation of each individual system, and 
= ^ytTp^'^^JtTpAp is the standard deviation of A. 

Now determine the amount of disturbance induced on the systems in the set. The 
state of any one of the systems in the set after the coupling with the pointer is given by 
tracing over all the other systems and the pointer state. Since the systems are identical 

by symmetry, look just at the first: 



p{t)=Y^ J dqdq''il;{q)i;{q')paa'\a){a'\p, 



0202 • ■ • raNO-N 



aa' a2---aN 

{q' + (a' + a2 + . . . + aN)lt\q + (a + 02 + . . . + aN)lt). (5) 



That is, the off-diagonal parts of pi are reduced by an amount 1 — / dqil^[q)il^{q-\-{a — a')'^f). 
A convenient initial pointer state IV') is a Gaussian wave packet centered at with standard 
deviation AQ (analogous to a coherent state of the electromagnetic field). In this case, it is 
easily seen that the effect of the coupling to the pointer is to multiply the aa! off-diagonal 
terms of p by a factor e-(T*(«-«'))V2AQ2_ rj,^^^ .^^ ^^^^ {^tAAf/2AQ^ e « 1, the 
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effect of coupling each member of the set to the same pointer is essentially the same as the 
effect of coupling each member to a different measuring apparatus, with a perturbation of 
size e = (7tA^/A(5)^. This model of measurement can be seen to be equivalent to the 
abstract POVM given above. 

It is interesting to note that the 'weakness' of this model of measurement can be 
tuned by adjusting the spread AQ of the initial pointer wave packet. As AQ becomes 
small, the measurement becomes stronger and stronger, revealing more information about 
an individual system while perturbing its wave function more and more. In the limit that 
AQ 0, this model reduces to von Neumann's original model of strong measurement. 
AQ acts as a knob that allows us to tune continuously from weak to strong measurement. 

We can weakly measure several observables simultaneously by adjoining several 
pointer variables and coupling ^^'jiOiQi. In the limit 7 — > 0, A?" ^ 00, the tth 
pointer provides an accurate assessment of (O^) while perturbing each system by as small 
an amount as desired. Note that need not commute with each other: in the weak 
measurement limit where the AQi are large and 7^ are small, the measurements do not 
interfere with eachother. By monitoring — 1 observables, one can obtain an assessment 
of all terms in the density matrix simultaneously. 

This concludes the detailed discussion of weak measurement. To summarize: by 
adjoining a suitable measuring apparatus and making the number of systems N in the set 
large, one can obtain the density matrix to a precision 5 = ^Jl/eN"^ + 1/N while perturbing 
each system by an amount e{dP — 1). By making N sufficiently large, the single-system 
reduced density matrix can be determined to arbitrary precision while perturbing each 
system by an arbitrarily small amount. It can be seen that the detailed model gives the 
same results as the abstract model of weak measurement given above. 

It is interesting to investigate whether the underlying statistics (Fermionic or Bosonic) 
of the systems in the set affect the results above. Since both wave functions and interac- 
tions are assumed to be symmetric, the results derived above hold equally well for Bosonic 
systems. If the systems are Fermionic, in contrast, they cannot be in completely identical 
states. However, if each system possesses additional degrees of freedom (position, for ex- 
ample, in the case of nuclear spins) that do not figure in the interaction with the measuring 
apparatus, then the discussion above applies to fermions as well. 

In fact, although we have assumed a symmetric situation in which the systems are 
described by identical density matrices p, this restriction is not necessary. If the systems 
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are prepared in the uncorrelated state pi <S> ■ ■ ■ <S> Pn, where in general pi ^ pj, then the 
entire set of results derived here applies to the determination and control of the average 
single-system density matrix p = (l/N) J2i Pi- 
Let us now assume that we can perform arbitrary weak measurements on a set of 
quantum systems, and feed the results of those measurements back continuously and co- 
herently using the well-known techniques of coherent control. That is, assume that we can 
implement arbitrary nonlinear unitary transformations as in equation (lb) and nonlinear 
Schrodinger equations as in equation (2). How might this technique be applied? 

The first potential use of this technique is simulation: a weak feedback controller could 
be used to simulate the dynamics of a variety of systems that obey a nonlinear Schrodinger 
equation. Nonlinear Schrodinger equations tend to arise in sets of weakly coupled quantum 
systems: as noted above in the context of the nonlinear Bloch equation, such coupled sets 
can be thought of as naturally occurring examples of weak feedback. For example, weak 
feedback can be used to simulate any set of systems that can be adequately described by 
a mean field theory, in which each system is coupled weakly to the expectation value of 
some operators on the ensemble as a whole. 

The use of a nonlinear Schrodinger equation is common in quantum optics to describe 
the evolution of photons that are weakly interacting with matter, as in a nonlinear optical 
fiber (23-24). As just noted, such an effect can be thought of as a naturally occurring ex- 
ample of quantum feedback by weak measurements: the atoms in the fiber weakly monitor 
and act on the photons. The use of weak feedback to create nonlinearities has the advan- 
tage that the form and strength of the nonlinearity induced by the quantum controller can 
be varied at will. For example, an optical weak feedback apparatus could be constructed 
by instrumenting a fiber with photodetectors and feeding back their signals to the fiber 
via electro-optic modulators (25). Such an optical controller could be used as a quantum 
analog computer to investigate the effect of time and spatially varying nonlinearities on 
the propagation of light down the fiber. It is important to note that such a fiber need not 
itself be nonlinear: all the nonlinearity could be supplied by the controller. In addition, 
weak feedback could be used to create and investigate the properties of optical solitons in 
a variety of nonlinear media. 

As noted by Haus (23), systems that obey nonlinear Schrodinger equations can be used 
to create Schrodinger's cats — quantum systems that exist in superpositions of two widely 
differing quasi-classical states. Although optical fibers are lossy and tend to introduce 
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decoherence, single-mode optical cavities of the sort constructed by Kimble et al. are 
good candidates for control by weak feedback (26). The very high Q of such cavities 
implies that the mode in the cavity is only weakly coupled to modes outside the cavity. 
Heterodyne monitoring of the cavity field therefore constitutes a weak measurement on 
the photons in the cavity. Nonlinearities induced by weak measurement could be used 
to create Schrodinger's cat states a la Haus. Wiseman and Milburn have also proposed a 
cavity quantum electrodynamics enactment of feedback via weak measurements to perform 
optical squeezing (7). 

Finally, quantum feedback via weak measurement could be used to create a novel form 
of quantum chaos. The usual linear Schrodinger equation does not exhibit sensitive depen- 
dence to initial conditions: the 'distance' between any two states |0) and as measured 
by their inner product remains constant (27-28). (Traditionally quantum chaos 

is not the study of sensitive dependence of quantum trajectories on initial conditions, 
but rather the study of quantized versions of classical chaotic systems.) The nonlinear 
Schrodinger equation, in contrast, need not preserve distances between quantum states, 
and can exhibit sensitive dependence on initial conditions (29-30). Quantum feedback via 
weak measurement, because it can be used to effect arbitrary nonlinear Schrodinger equa- 
tions, offers unique opportunities for investigating the sensitive dependence of quantum 
trajectories on initial conditions. Such ensemble quantum chaos could be used for example 
to construct a 'Schrodinger microscope' to detect and amplify small differences in quantum 
wave functions. 

We close by examining the relationship between nonlinearity induced by feedback 
of weak measurements and intrinsically nonlinear quantum mechanics. (Once again, the 
nonlinearity discussed in this paper is an effective nonlinearity: the underlying quantum 
dynamics of weak feedback is linear.) Nonlinear Schrodinger equations of the type found in 
equation (2) above are common in nonlinear quantum mechanics (31). Nonlinear quantum 
mechanics is known to exhibit a number of pathologies, including superluminal commu- 
nication (32), violations of the second law of thermodynamics (33), and the ability to 
solve hard computational problems (34). Since nonlinearity induced by weak feedback 
occurs entirely within the conventional framework of quantum mechanics, it cannot ex- 
hibit the first two of these pathologies. It might allow the solution of hard computational 
problems by the mechanism of the previous paragraph, viz., constructing a 'Schrodinger 
microscope' to detect small perturbations in the wave function of a quantum computer 
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(see also (35)); however, to obtain an exponential speed-up over classical computations 
one needs to amplify exponentially small differences in the wave function which in turn 
requires an exponentially large number of systems in the set. Nonetheless, it may be the 
case that weak feedback can be used to provide polynomial speed-ups to computational 
problems. 
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